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In this paper, we consider a reduced supersymmetric Yang-Mills integral with four super- 
charges by using a Gaussian approximation scheme and its improved version. We calculate 



the free energy and the expectation values of Polyakov loop and Wilson loop operators by 



extending the method employed in the bosonic case in the previous paper. Our results nicely 
match to the exact and the numerical results obtained before. The loop amplitudes exhibit 
good scaling behaviors similarly as in the bosonic case. The 't Hooft like large N limit leads 
simple formulas for the case of the loop length smaller. Also, the Polyakov loop and the 
Wilson loop are computed for the case of the loop length sufficiently large, where we see that 
the behavior of the Wilson loop reproduces the result simulated for a few smaller values of 
N at least qualitatively. 
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1 Introduction 



In this paper, we consider a supersymmetric Yang-Mills integral, which is the complete 
dimensional reduction of the TV = 1 supersymmetric Yang-Mills theory on the Euclidean 
four-dimensions with the gauge group SU(N). The classical action is S = + S^ F \ 
whereQ 



S {B) = -^tr[X^X v } 2 , = -^tr^(IV) Q/3 [X^]. (1.1) 

The variables (/z = 1, • ■ ■ , 4) are NxN traceless hermitian matrices, and ip a , xp a (a = 1, 2) 
are Grassmann valued traceless NxN complex matrices. For the Gamma matrices we 
take a representation Tk = i&k, = I2, with (k = 1, 2, 3) being the Pauli matrices. The 
partition function is defined as 



(dX)(d^) e~ s (1.2) 



with the integration measures 



4 n 2 -i dX A 2 ' n2 



(dx) = nn^ = n n • ( l3 ) 

fi=l A=l v Zn a=l \A=1 j 

The variables are expanded by a basis T A normalized by tr (T A T B ) = \S AB . 

The extended version of this model with sixteen supercharges has been proposed as a 
nonperturbative definition of IIB superstring theory ||]. It has no expansion parameters 
because the coupling constant g in the action can be absorbed away into an appropriate 
rescaling of the variables. So it is hard to analyze it, as usual, in the form of an expansion 
of some parameters. Under this situation, numerical calculations have been performed vig- 
orously 0, |3], |], ||, |fj . Analytic exact results for the partition function and the correlator of 
some operators have been obtained by using the cohomological field theory approach based 
on supersymmetry 0, ||, |9|0. 

In this paper, we perform analytic calculations expanding the methods of Gaussian and 
mean field approximations, discussed in the bosonic case in ref . |T7| , into the supersymmetric 



case. In the bosonic case, the partition function and various correlators including Polyakov 
and Wilson loops have been computed. The results for the Gaussian approximation fit 



2 This fermionic action differs from the action defined in eq. (2.1) in ref. || by the factor i. 

Unfortunately, the proof of the positivity of the fermion determinant A4 in H has a minor error. The 
identity o^-Mo^ = M* does not hold because A afi defined in Appendix A in |j[ is not real. However, the 
outline of the proof is correct. In fact, starting with our action and repeating the argument presented there, 
we can prove the positivity. Also, in the N = 2 case we are able to check it by an explicit calculation of the 
fermion determinant, which gives the answer identical to eq. (16) in the first paper of ref. [pj. Eventually, 
since the authors in |?| calculate everything by using the hermitian matrix T> = Jvft M, the factor i does not 
contribute and their final results are correct. 

3 For analytical results from other methods, see |l^, |ll|, Q. Also, for the cohomological field theory 
approach to quantum mechanical and field theoretical matrix models, see |l^, [l4|, |l5l fL6| . 
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better to the known results obtained numerically or analytically as increasing the number 
of the variables (the rank of the gauge group N or the dimension D). In particular, it can 
be seen that the results for the Polyakov and Wilson loops nicely reproduce the numerical 
results in ref. [0] for the region of the loop length L smaller. Also, we tried evaluating the 
asymptotic behavior of the loop amplitudes in the limit L — > oo by improving the method. 
Here, we perform a similar analysis for the supersymmetric case. For various supersymmetric 
quantum mechanical systems, the method of the Gaussian approximation has been examined 
in refs. [[[8], pT9f| , which suggests that it is necessary to formulate the Gaussian action so as to 
respect supersymmetry. This treatment is crucial, in particular, for observing the dynamical 
supersymmetry breaking from the behavior of the free energy. This is a reason why the 
Gaussian action is considered in the superfield formalism including auxiliary fields in refs. 
fl8| , |19H . We follow this lesson and consider the supersymmetric Gaussian action. In our case 
(D — 4), however we can easily realize the supersymmetry without using auxiliary fields. 
The situation is special for the four dimensions, so we will have to introduce auxiliary fields 
in the case of D = 6 or 10 (with eight or sixteen supercharges) ||20|| . 

This paper is organized as follows. In the next section, we consider the Gaussian ap- 
proximation to the supersymmetric Yang-Mills integral by extending the method used in the 
bosonic case |T7| so as to respect supersymmetry. We obtain the results for the partition 
function and the amplitudes of Polyakov and Wilson loops, which are reliable for N larger 
and for the loop length L smaller. In section 3, in order to get some insight for behaviors of 
the loop amplitudes for L large, we consider an improved version of the Gaussian approx- 
imation. The last section is devoted to summarize our results and mention possible future 
directions. 



2 Gaussian Approximation 

We consider the Gaussian approximation for the supersymmetric Yang-Mills integral (|1 . 1|) . 
As discussed in refs. [jTSL B|, respecting supersymmetry is crucial to realize a reasonable 



approximation. We start with the following Gaussian action: S = Sq B ^ + Sq F \ 



S { B) = ^tr (X M ) 2 , sf> = ^e aP tr (^ - ^). (2.1) 



This is supersymmetric, and the bosonic degrees of freedom and the fermionic ones are 
balanced. In fact, it is easy to see that the expectation value (S ) under the classical action 
S vanishes. In refs. |1| |19| , for the Gaussian approximation of supersymmetric quantum 
mechanical systems it is emphasized introducing auxiliary fields in order to balance the 
bosonic and fermionic degrees of freedom in Gaussian actions. It is equivalent to matching 
the freedom in the off-shell level because in the Gaussian action auxiliary fields contribute 
to the freedom as well as the dynamical variables. This is a reason why the Gaussian action 
is considered starting with the superfield formalism in refs. ]TB|, |19j. It is noted that in 
our zero-dimensional case the bosonic and fermionic degrees of freedom are already matched 
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without introducing auxiliary fields. In this approximation, the partition function (the free 
energy F) and the correlator of the operator A are expanded as 

F = F -J2 [ -^((S-So) n }c,o, (2.2) 

— Ti 

71=1 

OO C_1^7l 

(^) = (A) + ^^-f-((S-S ) n A) Cl o, (2-3) 

71=1 

where Fq and (• • -)o denote the free energy and the expectation value in the Gaussian theory. 
The suffix C means the connected part of the correlator. 



2.1 Partition Function 



First we consider the partition function in the Gaussian approximation. The outline is same 
as in the bosonic case ]T7| , but in the supersymmetric case to get a reasonable result we have 



to reorganize the series of S—So in fl2.2p so that the expansion respects the supersymmetry. In 
order to do so, we introduce a number we will call "order" : 1 assig ned to SW, S { B) and4 F) , 
and 1/2 to S^. It is easy to see from a diagramatic consideration as for the gap equations in 
ref. |18| that the expansion with respect to the total number of the order is natural from the 
point of view of the supersymmetry of the Gaussian action Sq. Straightforward calculations 
lead 



a 



(S (B) - S ) 
-i((^)-S ) 2 )c,o 



U(S {F) ) 4 )c,o 



4! 



-2(A^-l)ln — 



3^ 



N(N 2 - 1) 



-- 4 N(N 2 - l)a 2 r 4 



9 

45a 8 



lQg 4 
3a 6 r 4 



q 4 

N 2 (N 2 - 1) + ——N(N 2 - 1) 



2c/ 6 
5cr 4 r 8 



N 2 (N 2 - 1) 
N 2 (N 2 - 1), 



9 

3oV 
9 A 



N{N 2 - 1) 



(2.4) 



where the order 1 contributions are the second and third lines, and the order 2 contributions 
the forth and lowers. Considering the gap equations from the contribution up to the order 
1, i.e. 



d 
d 



F +(SW-S ) -~((SW) 2 ) C>0 
F + (SW-So)o-l((S^f)c,o 



0. 

= 0, 



(2.5) 
(2.6) 
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the parameters are determined as 

N / 

Then the Gaussian approximation for the free energy up to the order 2 is 




- 2 = \ 9 ^ r>= 9 -\ . (2.7) 



F= 1 -(N 2 -1) 



ln(Ng 2 ) + 1 h (corrections of (order) > 3) 



(2.8) 



where the second and third terms represent the corrections of the order 1 and 2, respectively. 
Let us compare this result with the exact one 0, 0] . The exact result is given by 

i,o, 1 9 Ar ( Ar+1 )/ 2 7r( Ar - 1 )/ 2 

a— = (V)-'""- 1 *^. r»= ■ (2 ' 9) 



and the large- N expanded form of the free energy becomes 

In M _ 1 

(2.10) 



-Fexact — — lnZ cxact — -j{N 2 — 1) 



, /»r 2\ 3 35 In N ^ . Ar 2 . 



Taking the ratio A = (F — F exact ) / F exa , ct for various values of N, we have A(N = 3) = 
-1.10408, A(N = 4) = -2.19021, A(N = 5) = -3.8156, A(N = 6) = -8.05504, A(N = 
100) = 0.464781, A(N = 1000) = 0.23609. Here we put g 2 = 1/2 for the same setting as 
in the bosonic case (eq. (4.6) in ref. flTT|). Note that differently from the bosonic case the 
inequality F < F + (S — S ) no longer holds because the integrals of Grassmann variables 
do not preserve the inequality. Thus A can take a negative value. Although for the smaller 
values of N (N = 3,4,5,6) the approximation does not make sense, for the larger values 
(N = 100, 1000) it becomes better as N increases. Also, for N larger it can be seen by 
comparing (|2.8|) with the expansion ( p.lOj ). On the other hand, in the bosonic case (see eq. 
(4.6) in ref. [|I7|) the approximation gives a reasonable result already at the smaller N. We 
can say that in the supersymmetric case the convergence of the Gaussian approximation for 
the partition function with respect to iV is slower than that in the bosonic case. 



2.2 Polyakov Loop 

Next we consider the expectation value of the Polyakov loop operator P(L) = -jjtr(e lLXl ). 
In the expansion ( |2.3|) , calculation of the correlators (P(L)) , —(P(L)(S — S ))c\o an d 
t;(P(L)(S( b > — So) 2 )c,o reduces to that of the purely bosonic contributions, and leads to 
the same formulas as obtained in the bosonic case |l7j], due to decoupling of the fermionic 
contributions or to the parity symmetry — > — of the Gaussian action. So our task 
up to the order 2 reduces essentially to evaluating the three correlators ^(P(L)(S^) 2 )c,o, 
- \ (P(L) (S^ - S ){S^) 2 )c,o and J (P(L) (S (F) ) 4 )c,o- In calculating the fermionic integrals, 
the formula 

(B(X)(S^) 2 ) C , = ^N{B{X)Si B) ) Cfi (2.11) 
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is useful. Here B(X) is a bosonic operator containing no ip a and ip a . The results are 
summarized as 



The lowest order 



(P(i)) = e -irMF 1-JV,2; 



-a 



(2.12) 



The order 1 correction 



■(P(L)(S^ - S )) Cfi + i(P(L)(^) 2 ) c ,o 



2<? 2 



-iV 



a 2 r 4 



LV 2 



2<? 4 



-N- 1 



e"^" 1 ) |(iV + 1)F (l - TV, 3; ^) 



-F 1-/V,2; 



(2.13) 



The order 2 correction 



liPimsW _ 5o)2)Ci0 _ ±<P(L)(S<*> - S )(S^f) c , + ^(F(^)(5 (F) ) 4 )c,o 



= e 4JV 



■(JV-l) 



3V 
2^ 2 



-TV 



a 2 r 4 



2^ 



-7V-1 + 



TV -3 /3a 8 



9' 



x 



2^ ) 
L 2 a T 



LV 4 

32iV 



+ e 4JV 



■(JV-l) 



'3a 4 a 2 r 4 V TV + 3 /3a 8 a 4 r 8 

r/V - — -TV - 1 + 



2g> 



2g* 



x 



4 V ^ 

7V + 3 ' 



2g 8 J 

F(l-N,5; 



4„4 



LV 
32iV 

L 2 a 2 \ 



■(JV-l) 



3o* 



(X T \ L a 



2g 8 J 327V 
L 2 a 2 \ 



TV \ / L 2 a 2 " 
3 ) F (*- N * — 



+ e 



a 4 r 8 



-(N-l) 



TV 
2 



9^ 

A9 2 



N-l 



N + 2 



3a 8 



2„2 S 



F 1-/V,4; 



L 2 a 



16^ 



■7V(3/V-2) + 



a 2 r 4 



iV^7V-2 N 



2^2 



N{10N - 1) 



L 2 a 



+ e" 



16c/ 8 v ~' ~ 7 J 4/V 



\2g* 
'9a 4 



2<? 4 

\ 2 



L 4 CT 4 

167V 2 



2„2 S 



F 2 - N, 3; 



L 2 a 



N-l 



3a 8 
16«?< 



■7V(37V + 2) 



^n(Kn 

9 A 



A9 2 



5 




L 2 a r 



(2.14) 



Let 



where F is the confluent hypergeometric function: F(a,f3; z) = J2^=o ^X^'-'-tp+n-i) nP 
us consider the value at the solution ( J2.7| ). Then, the order 1 correction vanishes, which 
is the same situation as in the bosonic case. Also, the 't Hooft like limit N — > oo with 
Ng 2 = G 2 fixed leads to the simple formula expressed by the Bessel functions: 

1 5 r 

(P(L)) = — JJ2y/a) + — J^y/a) - 4J 2 (2v^) + (corrections of (order) > 3) (2.15) 
la 16 L J 



where a = L 2 G/2. In ref. || some numerical results about Polyakov and Wilson loops 
are reported. We can compare our result with the numerical one presented there. Let us 
take Ng 2 = 48 and consider the quantity for various values of N keeping Ng 2 fixed, which 
gives the same setting as in ref. ||. Figure 5 in ref. || is the result to be compared with 
oursf]. The variable k/-Jg~ in the horizontal axis there corresponds to L in our setting. We 
plot our results up to the order 2 for various values of N in Fig. 1. The result exhibits 
a quite nice scaling behavior in the region < L < 6 in particular for N > 48 so that 
the curves lie almost on top of each other. In Fig. 2 we show that the approximation gets 
better as increasing the order from 1 to 2. Comparing the curves for the order 1 (the dotted 
line) and the order 2 (the solid line) with the numerical result (the gray dots), we see that 
the order 2 result nicely reproduces the numerical result for the region of L smaller (up to 
about 0.6)0 For the region of L larger the contributions of the higher orders are considered 
to become more important. For the order 2 result in the region of L larger, the amplitude 
of the oscillation is larger compared to the bosonic case (the gray line). Agreement for the 
region seems to be not so better than that in the bosonic long as comparing the 

results up to the order 2. 



2.3 Wilson Loop 

We calculate the expectation value of the Wilson loop W(L) = -Ltr ( e iLx le iLX 2e -iLX le -iLx 2 ^ 
by the Gaussian approximation up to the order 1. The formulas for (W(L)) and — (W(L)(S— 

4 There is a typo in the horizontal axis in Figure 5 in ref. Q. The horizontal axis represents k/^/g not 
k 2 /g. I thank authors of ref. || for informing me of this fact. 

5 The exact result of the Polyakov loop amplitude is non-analytic at L = because its second derivative 
(-^■tr(Xi) 2 ) is logarithmically divergent fl2|| . The Gaussian approximation up to the order 2 can not repro- 
duce this non-analyticity, but the approximation is considered to get better in the 't Hooft limit from the 
observation that the non-analyticity vanishes in the 't Hooft limit |2J] . 
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Figure 1: Polyakov loops up to the order 2 in the case N = 
gray line shows the case N — 16, the dotted line N = 48, 
dot-dashed line N = 400, and the solid line N = oo. 



-- 16, 48, 100, 400 and oo. The 
the dashed line N = 100, the 




Figure 2: Polyakov loop in the case N = oo up to the order 1 (the dotted line) and up to 
the order 2 (the solid line). The gray dots indicate some points of the numerical data for the 
N = 48 case in Figure 5 in ref. ||. The gray line is the result of the Gaussian approximation 
(up to the second order) for the bosonic case (D, N) = (4, oo) in Figure 3 in ref. [[L7j . 
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Figure 3: Wilson loops up to the order 1 for N = 16, 48 and oo. The dotted line shows the 
case N = 16, the dashed line N = 48, and the solid line iV = oo. The scaling is so good that 
the curves for the N = 48 and oo cases are on top of each other. 



S ))c,o are same as the ones obtained in the bosonic case [l7j (eqs. from (4.28) to (4.37) in 
the paper). Thus the quantity we have to evaluate is only ^(W(L)(S^) 2 )c,o, the result for 
which is 
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N 2 - 1 



f 2 (L)(-l+u(L)). 



(2.16) 



We used the functions f2(L) and u(L) defined in ref. |T7| (See eqs. (B.10) and (B.ll) in the 
paper). The explicit formula is too lengthy to write down directly. However, in the 't Hooft 
like limit it becomes remarkably simplified as 

(W(L)) = l-(-(Ji (2v^)) 2 -l 

+—. (9(Ji(2 v / o : )) 2 (J3(2v / ^)) 2 + 8(J 2 (2v^)) 4 ) + (corrections of (order) > 2)(2.17) 

A nice 
Our result can 



at" 

In Fig. 3, we plot the behavior for various values of iV with Ng 2 = 48 being fixed 
scaling behavior for N > 48 can be seen as well as in the bosonic case 
be compared with the numerical result in Figure 3 in ref. |§. Fig. 4 shows our result with 
some data for iV = 48 case in Figure 3 in ref. H . Taking into account the order 1 correction, 
the fitting of our result to the numerical one becomes better for the region of L smaller (up 
to about 0.7). Also, for the Wilson loop the behavior (the solid line) does not exhibit major 
discrepancies from that in the bosonic case (the gray line). 



3 Improved Mean Field Approximation 

In order to obtain some insight for the behavior of the Polyakov and Wilson loops in the 
region of L large, we consider about the improved mean field approximation along the same 



S 



Figure 4: Wilson loop in the case N = oo up to the leading order (the dotted line) and up 
to the order 1 (the solid line). The gray dots indicate some points of the numerical result 
for N = 48 in Figure 3 in ref. ||. The gray line, the result of the Gaussian approximation 
(up to the first order) for the bosonic case (D, N) = (4, oo) in Figure 6 in ref. [H], is added 
for comparison. 



line as discussed in the bosonic case | 17| . First, we begin with obtaining an interpretation 
for the Gaussian approximation as a mean field approximation. We take the following mean 
field action: Sm = Sffl + Sffi + X> where 



>M 



1 



q(F) 



— - — — tr(X u ) 2 ) 

Ag 2 N 2 -TJV V >* J 1 



\e a p\x {ipatpp) - -Leaptxtyaipp), 



+ 



N 



2g* (N 2 -1)(N 2 -2, 
1 1 1 
T 2 ~ 2g^(N 2 -1){N 2 -2 
1 1 1 



f 2 



2g A (N 2 - l)(N 2 -2) 



(tr('0«'0 /9)tr (ip-ylps)) + (tT(l/) a ^)tT(i)p7/) S )) , 

(tr(X M ) 2 tr(^/3)) - (tr(X^ a )tr(X^ /3 ))] , 
(tr(^) 2 tr(^)) " (tr(X^ Q )tr(X^))l , (3.1) 



A is a constant determined later. The bosonic mean field action comes from the con- 
traction of two bosons in and of four fermions in -(S^) 2 . The fermionic action Sm 
from the contraction of two bosons and two fermions in ^(S^) 2 . From the 0(4)-rotational 
symmetry and the SU (iV)-gauge symmetry, the following index structures are assumed: 



\{Xp)ij(Xv)kl{i)a)mnWp)pq) 



£ a fi£>yfi \Q'1'ijkl'-Ernnpq ^-^ijmn-^klpq b1ijpq1klmn\ j ( 3-3 ) 
fifiv^-aP ^ijkl^mnpq ~\~ b ' lijmn-^-klpq "I - b 1ijpql-klmn\ j (3.4) 
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where c, a, b, a' and b' are scalar constants with no indices, and Zy^ = SuSjk — jrSijSki- We 
evaluate (■ • •) in eqs. (|3~T| ) by replacing it with the expectation value under the mean field 
action (• • -) M , leaving the self-consistency conditions 

1 3N 2 N 4 1 N 2 2 
<j 2g z 2g^ t z g 4 

Here we take r to be real because we can always absorb the phases of r and f by rotating 
the fermions in Sm as ip — > e l6 tp, ip — > e~ ie tp. Eqs. (|3.5|) lead the solution (|2.7|). As same 



in the bosonic case |17), A is determined by the requirement Z = Zm (Zm is the partition 
function for the mean field action Sm-) as 

^ = -Y. [ —r((s-s' M ) n )c,M, (3.6) 

A 71=1 U - 

with S' M = Sm + S]p . Then, the partition function becomes 

2 oo (_-i)n 

F = - In Z = - In Z M = -2(N 2 - 1) In — - £ - ^)") CjM , (3.7) 

that reproduces the Gaussian approximation result in section 2.1. 

Next, we consider the unnormalized correlator (O)' = J(dX)(dipdiJj) Oe~ s where O is 
isotropic in the four dimensions and S77(iV)-invariant. We repeat the above argument re- 
garding (O)' as the partition function of a theory with the Boltzmann weight Oe~ s . Then, 
for such O, the correlators 

({X^)ij{X v )kiO), ((iJ a )ij(4>/3)kl(4> 1 )mn(4>s) Pq O) , {(X^ij(X v ) k l{^mntyp)pqO) 



have the index structures same as the amplitudes without O in ( |3.2| ), ( pT3| ) and (|3.4|), re- 
spectively. When starting with the mean field action: Sm = S M + S M + i, 

^M^ = ^ tr (^) 2 5 &aP = ^aP^a^Pp ~ ^£ a ptr4>a4>p, (3.8) 
(J 7" 

the self-consistency conditions give the formulas for a, f and f which are in the form same 
as eqs. (|3.1|) with the every correlators (• • •) replaced to (• • -) M = (■ ■ ■ 0)m/{0)m- In 
particular, in the case that O depends only on X^s, which is a situation we will consider 
from now, these conditions read 

- - 3 ^ (\r(X?)° N2 ~ 1 1 

a 2 ^A^-lV^ ) M N (i tr(X/i)2) o > 

h = wi^i^ tIix ^ )Z - (3 - 9) 

Note that the computation of the quantity (-^:tr(X M ) 2 )^ / completely reduces to the one in 
the bosonic case, because the fermionic contribution of the mean field action is factorized 
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and cancelled in each expectation value (• • -)m- From the requirement (O)' = (0)' M , after 
the same calculation as in the bosonic case [17], A is determined as 



1 

A 



-E 



m=l 



ml 



((s-s' m Dc,m-^ 



1 + 



E ^zt-(o(s - s' M ) n ) 



CM 



(3.10) 



The expression of the correlator becomes 

(oy- 



„ 2 x 2(JV 2 -1) 



-1/A 



f 2 



(O) 



(3-11) 



Finally, dividing this by the partition function evaluated up to the order same as done for 
(O)', the formula of the correlator (O) is obtained. In the following two subsections, we try 
to evaluate the amplitudes of the Polyakov and Wilson loops for the case of L large. 



3.1 Polyakov Loop 



For the Polyakov loop amplitude, we take the operator O as O 



iLX u 



1 1 

2-4 N 



Ej=i tr(e lLX » + 



Obviously, if the rotational symmetry of 0(4) is not sponateously brokenP], (O) is 



equal to (P(L)). When x = L 2 a 2 and y = L 3 r 2 large, the consistency conditions 
be solved in the expanded form: 



can 



x± 

y± 



, ± --' : < ( -j- (N + l)L 4 g 2 ) + ^N(N + 2) + O^g 2 )- 1 / 3 ), 



e ±W6 (18(iV+1)) l/6 (L Y ) 2/3 
1/2 



1 



1 + -NUN + 5)e Ti " /3 ( — — — 

3 V ; \16N+lL 4 g 2 



1/3 



+ 0((L 4 ^ 



4 2\-2/3 



(3.12) 



We do not consider the solution such that Rex < 0, because it leads to the unphysical result 
blowing up as L — >• oo. Corresponding to the above x, A is given by 



if (N~l)(^x + N + 0(x- 
A I (JV- 1) (^x + 3N + 2- ^(N + l) 2 + 0(x 



-i 



(up to the order 1) 
(up to the order 2). 



(3.13) 



If we regard the solution as the sum of the ± branches similarly to the bosonic case [JL7|] 



dividing the result corresponding to ( 3.11 ) by the partition function up to the order 1 or 2, 



6 At least, in the case of N finite, this always holds. When it is spontaneously broken, if we regard the 
correlator (• • •) as the sum of the expectation value at each superselection sector, the equality (O) = (P(L)) 
holds. The spontaneous breaking of the rotational (Lorentz) symmetry in reduced Yang-Mills integrals is 
discussed in refs. g, ^ g || §| |§ . 
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we get the final answer 

(P(L)) ~ C H (LY)- {N - 1)2/6 



X cos 



exp 
N- 1 /16 



— w 



(<<) 



(n)" 



iV 



1/3 



7T 



3 (iV + l)LVJ --(iV-l)(iV + 2] 



where up to the order 1 (n = 1): = 1/4, and up to the order 2 (n = 2): U( 2 ) 
overall constants are 



, (3.14) 
3/8. The 



a 



(i) 



a 



(2) 



(-l^N-l N (N 2 -l)/2 

2^-2^! 
2^-^! 



(jV + 1 )(JV-l)(JV+2)/3 



12 g^(iV 2 -l)/3 ^ 16 x (AT-l)/3 



-(AT-l)(77V+5)/6 



17 V 3 / 

(iV + l)(JV-l)(iV + 2)/3 ^y^^ 73 f^y" _1)/ %(^-l)[&(tf+l)-?5] 



3 -; 



It is interesting to compare this with the exact answer obtained in ref. In |Tl] the 

eigenvalue distribution p(A) = (-^tr<5(A — Xi)) is exactly evaluated for the N = 2 case 
(see eq. (12) in ref. fill). The Polyakov loop is related to this quantity via the Fourier 
transformation: (P(L)) = J^ c oc dXe' lLX p(X). This integral can be easily estimated in the case 
L large, and we end up with 



(P(L)) 



r(§) \ 2 7 



exp 



2 I 2 7 J 



x [l + 0((L 4 s 



(3.15) 



Although our result reproduces the power (L 4 g 2 ) 1//3 of the exponential decay, the oscillating 
behavior can not be seen in the exact result, also the power behavior is different between 



the improved mean field result ((L 4 g 2 ) 1 / 6 ) and the exact result ((L 4 g 2 ] 
exponential decay behavior of the improved mean field result goes as exp 



1 / 3 ). Moreover, the 
-2 2 / 3 (2- 7 LV) 1/3 



for up to the order 1 and exp 



for up to the order 2. Since the co- 



-3 ■ 2- 1 / 3 (2- 7 LV)V3 

efficients are 2 2 / 3 = 1.58 • ■ ■ and 3 ■ 2 -1 / 3 = 2.38 ■ ■ -, we can not say that the approximation 
tends to converge to the exact results as increasing the order. About this disagreement, we 
could consider as follows. The Gaussian or mean field approximation becomes better when 
the number of variables is larger, which means that iV is larger. As is seen in the analysis of 
the partition function in section 2.1, the result begins exhibiting the nice convergence for iV 
larger than about 100. For smaller iV (N = 0(1)), the approximation does not make sense. 
There will be a possibility that the Polyakov loop has the same tendency as the partition 
function for the dependence of N and that our result makes sense for N large. In the result 
of the Gaussian approximation in section 2.2, because all the confluent hypergeometric func- 
tions appearing there are polynomials of L 2 <r 2 for N > 4, the oscillating behavior terminates 
at some finite value of L. When N increases, the degree of the polynomials increases and 
the oscillating behavior continues to larger value of L. Thus it seems to be not ridiculous 
to expect that the asymptotic behavior of the Polyakov loop has the oscillating part in the 
case of N large. 



12 



3.2 Wilson Loop 

We evaluate the Wilson loop amplitude up to the order 1. For the operator 

O = — Y —tT(e iLX ^e iLX "e- iLX ^e~ iLX "), 
4-3 V N 

we apply the procedure at the beginning of this section. The solution of the consistency 
condition (fOl) is given by 



a 



1 



f 2 = r 2 U 



3(JV + 1)!(JV-1)! 
1 1 



e 2 



3(iV + l)!(JV-l)!" \ 2 

1/4 



/ r2/T 2\2AT-2 ^ 

f^-j +0(e-^ 2 (LVr- 4 ) 
0(e- LV2 (L 2 a 



^^2^-2 

e 2 



2\4V-4 



.(,3.16) 



with a 2 = and r 2 = (jjpj . In the limit L —>■ oo, cr and r go to the values obtained by 
the Gaussian approximation. A is determined as 



i = V-D i - 



2 f _L_ 

iV!(JV-2)!\ 1 + iV 2 



2 | 2 (iV + l)(iV + 4) i 



3(iV + 3) (iV + 2)(iV + 3) 



[ e -^<T 2^2^-3' 



(lV) 2JV - 3 V 



After all, we obtain the final result up to the order 1 



2 a / 7-2^2X2^-2 



e 2 



(3.17) 



(W(L)> 



iV + 1 



exp 



N 2 - 1 
2iV 2 (jV + 3) 



1 2 [ 2 (iV + l)(iV + 4) ^ 2 



2 2 /r2^2\ 2iV-2 
L J g- J [ Li O \ 



N\(N-2)\N 2 { (N + 2)(N + 3) 



3(iV + 3) 



e 2 



O e-^(LV 



2 2\27V-3 



(3.18) 



This behavior has the same tendency as in the bosonic case [|17]]. Namely, (W(L)) has a 



nonzero 

2 

AT+1 



exp 



imit as L — > oo, and approaches to this limit from the below. The value of the limit 



2N 2 (N+3) 



becomes tiny as increasing N. It is consistent to the numerical result for 
L larger (up to L — 40) and N smaller (N = 2, 4, 8) in Figure 3(a) in ref. [O]. Since our 
approximation is expected to becomes better for N larger, we could say that for the Wilson 
loop our result nicely reproduces the true behavior at least qualitatively. 
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4 Discussions 



We have considered the reduced four-dimensional supersymmetric Yang-Mills intergral with 
the four supercharges by using the Gaussian and mean field approximations. We extended 
the approximation scheme considered in the bosonic Yang-Mills integrals so that the super- 
symmetry is respected. As a result, we obtained reasonable results for the partition function 
and the correlators of the Polyakov and Wilson loops. The result of the partition function 
shows a nice convergence into the exact value as N goes larger. Also the Gaussian approx- 
imation for the Polyakov and Wilson loops fits well the numerical results for the region of 
the loop length L smaller. Next we considered the improved version of the mean field ap- 
proximation in order to analyze the behavior of the loop amplitudes for L large. Although 
the result of the Polyakov loop does not exhibit good agreement to the exact result for the 
N = 2 case, it could be expected that the situation gets better for N large as we saw in the 
analysis of the partition function. The asymptotic behavior of the Wilson loop shows a nice 
behavior that reproduces the numerical result at least qualitatively. 

Comparing with the situation in the bosonic case, we notice these two things. For both of 
bosonic and supersymmetric N goes larger, the result of the partition function tends 

to converge into the region where the approximation is reliable. For the supersymmetric case, 
however it is slower than that in the bosonic case. Also, in the Gaussian approximation for 
the Polyakov loop, the difference for the region L > 1 between the approximate result and 
the numerical one is larger than in the bosonic case. Thus, we can conclude that as N goes 
smaller and L larger the Gaussian approximation gets more worse comparing with the result 
in the bosonic case as long as looking at the first few orders. Since we have examined the first 
few terms in the expansion around the Gaussian action, it is obvious to need investigate the 
higher order corrections. At present, however, computation up to the higher order gets much 
more complicated than that for the first few orders. We would need a calculation technique 
making the computation greatly simplified or an improvement of the approximation so as to 
give better results up to the first few orders. 

An interesting direction is an investigation of the case of the supersymmetric Yang-Mills 
integral with eight or sixteen supercharges. In this case, there is a technical difficulty for 
obtaining reliable results by numerical simulations due to the fermion determinant being not 
positive definite^. Since our method seems to be extended also for this case by introducing 
auxiliary fields, it would be expected to give concrete results 
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7 There are some attempts to find a way out of this problem in refs. 
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